On Stable Pareto Laws 
in a Hierarchical Model of Economy* 

A. M. Chebotarev 

Quantum Statistics Department, Moscow State University^ 
Moscow 119899, Russia 

Abstract 

This study considers a model of the income distribution of agents whose pairwise 
interaction is asymmetric and price- invariant. Asymmetric transactions are typical 
for chain-trading groups who arrange their business such that commodities move 
from senior to junior partners and money moves in the opposite direction. The price- 
invariance of transactions means that the probability of a pairwise interaction is a 
function of the ratio of incomes, which is independent of the price scale or absolute 
income level. These two features characterize the hierarchical model. The income 
distribution in this class of models is a well-defined double-Pareto function, which 
possesses Pareto tails for the upper and lower incomes. For gross and net upper 
incomes, the model predicts definite values of the Pareto exponents, Cgross and Onet, 
which are stable with respect to quantitative variation of the pair-interaction. The 
Pareto exponents are also stable with respect to the choice of a demand function 
within two classes of status-dependent behavior of agents: linear demand (ogross = 
1, Cnet = 2) and unlimited slowly varying demand (agross = «net = !)• For the 
sigmoidal demand that describes limited returns, Ogross = «nct = 1 + a, with some 
a > satisfying a transcendental equation. The low-income distribution may be 
singular or vanishing in the neighborhood of the minimal income; in any case, it is 
-Li-integrable and its Pareto exponent is given explicitly. 

The theory used in the present study is based on a simple balance equation and 
new results from multiplicative Markov chains and exponential moments of random 
geometric progressions. 



1 Introduction 

1.1 A brief historical survey 

The orthodox interpretation of statistical distributions in economics assumes a certain 
stochastic process whose invariant probability measure characterizes the given distribu- 
tion. Champernowne [1] made a remarkable attempt to model the Paretian distribu- 
tion, and obtained several discrete versions of the distribution as stationary probability 

*This work is supported by DFG 436 RUS 113/779/0-1. 
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measures for a class of Markov chains on half-line. Starting from an empirical analysis 
of UK tax records for the period 1951-52, Champernowne reconstructed a (12 x 15)- 
approximation of the stochastic matrix of the Markov chain, and noted a definite "ten- 
dency for the lowest income to shift upwards" [1]. In a sense, his basic assumptions 
(asymmetry and homogeneity of the transition probability) are inherited by the model 
described in the present study. To ensure the probabilistic interpretation and to obtain 
a closed equation for the stationary distribution, Champernowne introduced a boundary 
condition for the occupation number of the lowest income range. This fare for mathemat- 
ical correctness was not reasonably explained in economic terms, but a similar model of 
multiplicative stochastic process with "reset" events (i.e., with stochastic renewal of the 
process at the boundary) [2] was successfully applied in [3] for a Monte Carlo simulation of 
variations in income distribution within Japan. In the present paper, only the continuity 
assumption is imposed at the boundary (at infinity and at the point corresponding to the 
minimal income) that defines a stationary density distribution. 

To fix notation, recall that the cumulative probability distribution Pt{x > s) of a 
random variable x has the Pareto tail with exponent a if Pr(x > s) = 0(s~"), a > 0, for 
large s^. 

A number of empirical facts in economics concerning the Pareto distribution can be 
qualitatively explained using generalizations of the Yule stochastic process [4], which 
was proposed in 1937 as a simple dynamical model of biological taxa. In view of a 
certain similarity between the present model and Yule's model, it is worth explaining the 
basics of Yule's model in a macroeconomic context. Suppose that the starting capital 
equals 1 and a nonrandom consumer demand provides the business community with fixed 
1 + m units of income per unit of time. Each step of the time-evolution consists of two 
events: (1) a new agent obtains the starting capital and joins the business community; 
(2) the remaining m units of income are randomly distributed among the existing agents 
according to their status-demand, for example, with partial probabilities proportional to 
the existing distribution of capital (linear demand). The construction of a stationary 
solution in a biological context is discussed in [5] . The cumulative stationary distribution 
has the power tail with the following exponent: 

1 

a — 1 + a, a — — . 

m 

If the business society only rarely recruits new agents, a is small and a ~ 1^. 

In addition to linear demand (^-demand), the class of slowly varying functions is 
considered. From an empirical viewpoint, in the linear case the size of a business is far 
from market limitations or the constraints of resources. The slowly varying demand (it- 
demand) carries features of utility and sigmoidal functions that describe the diminishing 
returns of investments and the effects of market or production constraints. 

^In fact, the conventional terminology is not completely correct: Zipf constructed his plots with the 
probability axis y, while Pareto used x for the same purpose [5]. 

^Let n be the number of time units per year. Then, n m w n (1 + m) = cr is the size of the market and 
1/ = n is the intensity of the creation of vacancies. If a > 1, it is considered that a = a — 1 is the horizon 
of the power tail, as the heaviest tail with a = 1 normally overshadows the remaining components of a 
sampling. In these terms. Yule's law states: "Horizon x market size ~ intensity of vacancies ' creation" . 
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Empirical studies of the top income distribution show that the Pareto exponent a of 
the cumulative distribution normally belongs in the interval [1, 3]. There is evidence for 
both stability [6]-[10] and instability [11, 12] for this indicator over time. In the latter two 
reviews, the authors discuss the relationship between tax legislation and the distribution 
of personal income. The reviews present variation of the Pareto exponent in the range 
[1.3, 3.5], based on tax records from Great Britain and the Netherlands for the years 
1907-1999, and records from the US, Canada, and France over shorter periods. The 
detailed records enable us to monitor details of income distribution for the top 0.05% to 
10% of earners. The Pareto index a = 1 for the net income of firms and corporations 
has been observed by many authors (see [9, 17, 18]), including the author of the present 
paper [19, 20]. These tails are the least informative in the sense that they overshadow 
the possibly important details of less heavy tails that contribute to the general picture. 




Figure 1: The two lower curves in the left-hand figure display the cumulative probabilities 
defined by 2000 Pareto-distributed random points with exponents a = 1 and a = 2, respectively, 
while the upper curve shows their 1:1 mixture: a^ix ~ 2.3. This illustrates the shadowing effect 
created by the heaviest tail. The point-plot on the right shows amix(a^) as a function of the share 
X of the heaviest tail varying from 5% to 95% over the total sampling of 4000 random points. 

A solvable model for the Pareto income distribution was recently proposed by Reed 
[13]. In fact. Reed's approach is based on the Poincare ergodic conjecture, which states 
that the time-average coincides with the average over the ensemble. The ensemble con- 
sidered by both Reed and Gabaix [14] consists of noninteracting agents whose income 
is described by geometric Brownian motion and whose life-time is an exponentially dis- 
tributed random variable. By averaging the individual income distribution over the ran- 
dom life-time, the double- Pareto distribution [13] is derived^. According to the Poincare 
conjecture, the obtained distribution describes the income distribution over the ensemble 
of ideal (i.e., noninteracting) agents. 

The present paper considers an ensemble of asymmetrically interacting agents. The 
balance equation that describes the stationary distribution of income is derived using 
arguments that are similar to those used by Bouchaud [15] for symmetric pair-interactions. 
The present equation enables the calculation of the pair of Pareto exponents a and d > —1 

correspondent informed me that in 1949, E. Fermi described a similar mechanism in an astrophysical 
context, but I have been unable to confirm this. 
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for the top and bottom income tails. The value of the Pareto exponent a^ct = 1? derived 
below for hierarchical models with w-demand, is normal for the top revenues of firms 
and corporations in Japan (see [3, 6]) and is in agreement with estimates for Europe, 
Britain, and the USA (see [10, 12, 13, 16]). Indirect reconstruction of the gross income 
distribution for the top 5% of householders, based on sale statistics for the national 
car market and direct estimates of the capital distribution among the 200 largest banks 
in Russia, yields the same value [19, 20]. A possible qualitative explanation for this 
rather general phenomena is the existence of a dominating corporate subculture of profit- 
making communities. The aim of this paper is to consider the extreme case of a corporate 
community: the hierarchical structure. A sketch of this approach was published in [21]. 

1.2 Assumptions and results 

In the business world, the hierarchical environment maintains and permanently repro- 
duces the specific circulation of money and commodities; corporations, network-marketing 
and chain-trading structures, shadow economy groups, etc., aim to organize supply and 
demand such that commodities and money move in opposite directions. Within the hi- 
erarchy, an agent consumes the commodities distributed by senior partners, while at the 
same time, he distributes them to junior agents [22]. 

Clearly, the above picture is not strictly true for agents with low income who are 
receiving money (e.g., low wages or pensions, welfare subsidies or grants) from the state. 
To simplify mathematical considerations, taxation and the redistribution policy of the 
state is ignored, and a toy model is considered that describes the stable distribution of 
unilateral money flow in a hierarchical structure under the assumption that the intensity 
TZ{s, x) of transactions between the pair of agents is an asymmetric and price-invariant 
function of incomes s and x of these agents; that is, Tl{s,x) = R{s/x) > if s > a;, and 
R{s/x) = if s < a;. These assumptions characterize this class of hierarchical models. 
We also assume that income is (mathematically) unlimited and that the minimal income 
So is positive, i.e., s e [sq, oo), Sq > 0. 

The money flow controlled by an agent with net income s consists of the welfare 
payment V{sq,s) (cash beneflts, social grants, etc.) and the contribution B{n,s) from 
subordinate agents. Here, n = n{s) e Li{so, oo) is the density of agents with net income 
s, and B{n, s) is a function of s and a functional of n. The hierarchy assumption means 
that B{n,s) depends on the values of n{x) restricted to x e [so,s], and B{n,s) — > 
as s — > Sq. The welfare 7^(50) s) decreases and vanishes for all sufficiently large s. The 
conjectured price-invariance of the model implies that 7^(501 s) = P{s/so). 

To maintain the social status, an agent with net income s pays out H = H{n, s) to 
senior partners. We assume that the hierarchy is free of competition; in other words, 
we conjecture that (at least in the neighborhood of the stable distribution) variations 
in the number of senior agents docs not affect consumer demand and price levels. This 
assumption implies that the individual consumption costs H — H[s) are independent of 
n. Thus, the net income is equal to 

s^P{s/so) + B{n,s)-H{s), (1) 

and the sum g{s) — P{s/so) + B{n, s) = s + H{s) is the gross income. Suppose that the 



4 



consumption costs H{s) decrease as s | Sq. Hence, the gross income g{s) — s + H[s) is 
an increasing function of s and the equation 

So: So + H{so) = P(l) (2) 

has a unique solution. Furthermore, (2) is regarded as an analytical definition of the 
minimal income Sq, while (1) is considered as an equation with respect to n(s). 

In Section 2, B and H are represented in terms of the pair- interaction intensity R{x) 
and the demand function a{s). The main statement of this section is a conditional 
proposition: if the solution n{s) of equation (1) can be represented in the form 

n{s)=x''{l-xYp{x)/a{x-^), x = so/s e {0,1] (3) 

where p is a continuous, strictly positive, and uniformly bounded function on [0, 1], then 

d=-l + ^>-l, C(l) = l-P'(l) + 2i?o>0, Ro''^ r R{x)dx. (4) 

G(ij Ji 

In addition, it is demonstrated that Ogross = 1, dnct = 2 for the linear demand, Ogross = 
anet = 1 for any slowly varying demand, and Ogross = c^net = 1 + a for the sigmoidal 
demand, where a is a unique solution of the transcendental equation 

a:-J^ ."i?(.)rf.= l + ^, ^o=lima(.). 

Section 3 considers the results of numerical simulation of p that ilhistrate various op- 
portunities of the model. The existence of the income distribution n(s) in the form (3) 
is proved in Section 4. p is first represented in terms of the pair of multiplicative Markov 
chains, and {rjn}, that describe the function p{s) either a.s s ^ oo or s ^ Sq. In 
Section 5, possible economical issues are discussed. Section 6 consists of a mathematical 
appendix that contains four basic lemmas on multiplicative Markov chains. It is antici- 
pated that mathematical facts considered in the last section will be useful for the readers 
who use multiplicative Markov processes in econophysics and mathematical finance. 



2 Balance equation and continuity conditions 

2.1 Balance equation 

Suppose that the money fiow from agents with income x to n{s) agents with higher income 
s is proportional to the density n{x), consumer demand a{x), and the pair-potential 
R{s/x) of interactions. As this income is shared among n(s) senior agents, the personal 
inflow functional B{n, s) can be written as follows: 

B(n, s) — . , / aix) Risl x) nix) dx. 
n{s) Jso 

The given rational agent behaves in exactly the same way toward the senior partners as 
his junior partners: the cost of his demand to agents with higher income x is proportional 
to the intensity of interactions and to his demand function a{s): 

p OO /"OO 

H(s)^a(s) J R(x/s)dx^sa(s)Ro, = R(y)dy. 
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Taking into account the welfare income P, the balance equation (1) is rewritten as follows: 
s — P{s/so)-\ — --— / a{x) R{s/ x)n{x) dx — a{s) / R{x/s)dx. (5) 

TtyS) J so Js 

This relation is considered as an equation for the stationary distribution n{s) of the net 
income s. The equation g = s {1 + (t{s) Rq) is an obvious statement on the relation 
between the gross and net income, g and s, which readily follows from (5). 

For further use of the model, it is convenient to rewrite (5) as a linear homogeneous 
integral equation: 

/ a{x) R{s/x) n{x) dx — (s — P{s) + s a{s)Ro) n{s). (6) 
J so ^ ' 

The choice of a smooth monotone demand a{s) < s is restricted to the following three 
classes: 

(i): a{s) — s for all s > 1, (s): cr(s) = s for s = 0(1); a{s) Sq as s — > oo, 

(u): a{s) = s for s = 0(1); 1, a{s) oo as s — > oo. 

The slowly varying functions a{s) = 1+lns, s > 1 (utility function), a{x) = Sq s/{So + s — 
1), and 5*0 > 1 (sigmoidal function) are examples of u- and s-demands. These functions 
describe various kinds of diminishing returns. 

Because the contribution from junior partners, B{n,s), vanishes as s J, sq, to ensure 
consistency in equations (5) and (6), equation (2) is rewritten as follows: 

so-P{l) + soa{so)Ro = Q (7) 

and consider the equation to be the definition of Sq- Equations (6) and (7) represent the 
analytical setup of the model. The aim is to describe the asymptotic behavior of n{s) 
as s — s> oo and s — > Sq and to prove the uniqueness of the solution. For simplicity of 
notation, the price scale is chosen such that sq = 1. 



2.2 Consequences of the continuity conjecture 

To study the problem on a compact set, the variable s G [1, oo) is changed to x' = 1/s G 
(0, 1] in (6) and (7). According to (7), the integral in the left-hand side of (6) vanishes at 
the point s — x — 1. Hence, the right-hand side of (6) can be rewritten as the following 
product: 

s - P(s) + s a(s) Ro = a{s) (s - 1) C(s-^) = x-^a{x~^)(l - x) C(x), (8) 

where a{x^^) = x~^ for x = 0(1). 

As s ^ oo and x = 1/s — > 0, the leading terms in the left- and right-hand sides 
of (8) arc x^^a{x^^)C{0) and x~^{l + a{x~^) Rq), respectively. Therefore, for C(0) = 
Rq + a{oo)~^, we have 

nfrw _ / -^0 for ^- s-nd li-demands, , , 

\ Ro + So^ fors-demand, 
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because cr(oo) = oo for £ and tt-demands, and cr(oo) = 5*0 for s-demand. As x —>■ 1, we 
obtain 

C(l) = 1 + 2i?o - > 1 + 2i?o > 0, (10) 

where P'{1) < because the welfare P{s) is normally a decreasing (or finite) function. 
Therefore, the left-hand side of (8) is an increasing function that vanishes only at s = 1. 
For this reason, it is supposed that c* > exists such that 

c{x) = C{x)/C{0) > Va; e (0, 1]. 

Suppose that P{x~^) vanishes for sufficiently small x, say for x e (0, x*], and a{s) — s for 
s e [1, l/x*] ^. It is clear that for x e (0,a;*], the assumption a{s) < s implies 

(1 - x)C{x) _ . 1 ^11 + x/Ro for i- and li-demands, , . 

Ro ~ (j{x-^)Ro ~ 1 1 + l/SoRo for s-demand. ^ 

Let us search for a positive solution of the balance equation (6) in the factorized form 

n{l/x)^x''{l-xYp{x)/a{x~'^), x e [0,1], b > 1, d > -1. (12) 

As the normalization of n is not essential, p(l) = 1 is used. Let us derive an equation for 
p and calculate a priori values of b and d. Set 

r(s) = i?(s)/i?o, m„ = ^°°sV(s)ds = 1, mo = 1 (13) 

and suppose that the probability density r{x) is a bounded continuous function, r(l) > 0. 

Theorem 1 The relation d — R{1)/C{1) — 1 is necessary for the continuity of a strictly 

positive bounded function p{x) at the point x = 1 for all continuous demand functions 
such that a(l) = 1. Continuity at the point a; = implies = 3, bu — 2, and bg > 2 for 
s-demand; more precisely, bg satisfies the transcendental equation 

bg-. / s'^-\{s)ds=l + ——. (14) 

Jl Do Ko 

Proof. According to equations (8) and (12), rewrite equation (6) in the variables x — 1/s, 
y = 1/^: 

n{s) (s - P(s) + s a{s) Ro) =Ro x^~\l - xy+^c{x) p{x) 

= r a{y)R{s/y)n{y)dy= [' a{z-^) z-^R{sz) n{z-') dz 

Jl Jl/s 

= j\^-\l-zYR{z/x)p{z)dz. 

Jx 

*A point a;* € (0, 1) will appear in some of the assumptions below. In fact, these additional assumptions 
can be imposed at different points, but a fixed point is used for simplicity of notations. 
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This integral equation for p[x) can be represented in two equivalent forms which are 
indepemdent of a: 

(15) 

1 6-2 /^l -a^l/"^ 



(1 - x)C{x) Ji 

If d > — 1 and X > is fixed, the measure //(x, dz) — {d+ 1)(1 — zY/{l — xY'^^ dz is a 
probability measure on {x, 1]. If p{x) is a bounded measurable function that is continuous 
at the point x = 1, so is the product 

p{x, z) = {z/xf~'^ R{z/x) p{z), z e {x, 1) 

and limp(a;, 2;) = limp(a;, ^) = as x — >• 1, 2; e {x,l]. Denote with the 

mathematical expectation with respect to the probability measure iJ,{x, dz). Hence, 

limE^p(x, •) — limEa;p(x, •) — R{1) 

Multiplying (15) by d + 1, we obtain the condition that is necessary for the continuity of 
P- 

(d+l)p{l) = \M^p{xr)/{xC{x)) = R{1) p{l)/C{l). (17) 

If p(l) ^ 0, this equation yields (4): d = i?(l)/C(l) - 1 = r(l)/c(l) - 1 > -1. 

As X ^ 0, the left-hand side of (16) clearly converges to p(0); the Lebesgue theorem 
on dominated convergence implies that the right-hand side of equation (16) converges to 



p(0) Ro 

C{0) 



mb-2 = p(0), rua =^ ^ s"r(s) ds. 



As C(0)/i?o = 1 for £- and w-demands, and C{0)/Ro = 1 + I/SqRq for s-demand, the 
necessary continuity condition at x = is the equation with respect to b 

{1 for i- and -u- demands, ^ . 

1 + 1/RoSo fors-demand ^ ^ 

provided p(0) 7^ (see (9) and (13)). The moment rua is an increasing function of a 
because r is a regular probability density on [1, 00); hence, equations (14) or (18), in each 
case, have a unique solution. □ 

Remark 1. As a{{xy)~^)/cr{x~^) ^ 1 as ,x — >■ 0, this result remains valid if p{x) = 
g > 0. This will be the case for w-demand. It is also clear that the Pareto exponent for 

n{s)0{cr{x~^)~^) is the same as that for n{s). 

Remark 2. li 6 = I/SqRo is a small constant and the measure r(s) ds has at least one 
moment, the equation nia =1-1-5 can be approximately solved: 



i.e., a = 1 + a, similar to the Yule case^. 

Thus, it is demonstrated that the continuity and strict positivity of n{s) implies the 
expansion: 

, _iN h/-, \rf / \ / / -i\ 7 -. f 1 for £- and M-demands, 

u{x )^x\l-xYp{x)/a{x ), i'-l^^net^l for .-demand. 

In addition, for all kinds of demands, equation (15) can finally be rewritten as follows: 

where a = for I- and M-demands, and C is the only function that is dependent on a 
(see equation (8)). 

2.3 Pareto exponents for gross income 

The relationship g = s + a{s) Rq between gross and net income allows one to calculate the 
Pareto exponents for the gross income distribution. First consider the ^-demand. Denote 
the unique positive root of the equation g — s + Rq with s{g), where g is gross income. 
The density distribution of g then equals Ne{g): 



Ne{g) = ne{s{g)) s\g), s{g) = l/^l + m,g = 0{g'/'), 

where ni{s) = 0{s~^), n^{s{g)) = 0{g~^^'^), and s'{g) = 0{g~^^'^). Finally, we obtain 
Ne{g) = 0{g~'^), i.e., the Pareto exponent of the gross income Ogross = Onet — 1 = 1. 

Now consider the it-demand. Because cr(s) is a slowly varying and increasing function, 
there exists a finite constant such that a{s) < a^s^ for any £ > 0. Now denote the 
unique positive root of the equation g = s {1 + (t{s) Rq) with s{g), and let N^ig) be the 
probability density of g. The monotonicity of cr(s) implies the following inequalities: 

\ 9 ^ 9 ^ 9 (0,0^ 

g-syg)- - i + - l + a,g-Ro ^ ' 

By definition of the cumulative probability, we have 

POO POO POO 

/ Nu{g)dg^ / nu{s{g)) s'{g) dg = / nu{g)dg. 

Jy Jy Js{y) 

As nu{g) = 0{g~^), the above inequality (20) yields the two-sided estimates for the 
cumulative probability fu{g > y}'- 

rco roo rco 

0{y-') = / n^{g) dg < / riuig) dg = / dg = ^"^ K{9 > y} 

Jy •'s[y) Jy 

roo 

< / riuig) dg = y~^0{l + a,y' Ro), 

Jy/{l+aey^ Ro) 



^Thc constant Rq characterizes an intensity of pairwisc interactions. In this context, I/-R0 is the 
average time between interactions. On the other hand, Sq can be interpreted as the volume of a "market" . 
The difference a — 1 between the Pareto exponent a (a > 1) and the exponent of the "heaviest" tail is 
termed the horizon of the distribution. In these terms, the Pareto exponent that corresponds to the 
sigmoidal demand satisfies the following law: '^Horizon x market volume ~ expectation- of-transaction 
time." 
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where Z is the normahzing constant and e is arbitrarily smaU. Hence, passing to the hmit 
as J/ — >• oo, £ — >• 0, we have 

Ogross = -lim ln(P{g' > y})/ Iny = -hm ln(P{g' > y})/ Iny = 1 = Onet- 

For s-demand, gross and net income have the same order: g((s) = s (1 + (t(s) Rq) = 
s + 0{s), as a{s) is bounded by Sq. Hence, 

1 f°° 

Onet = Ogross = 1 + 0:, « : 1 + -— - = / r{s) ds. (21) 

OqUq Jl 

The above estimates of b and d define meaningful asymptotical properties of n{s) 
provided that p{x) is bounded from above and below. The proof of this fact is a nontrivial 
mathematical problem that will be solved via a probabilistic representation of the solution 
of equations (6) and (7); however, first the results of an empirical study of this problem 
are considered. 

3 Numerical simulation of p 

With respect to the analytical proof of the boundedness and the strict positivity of p{x) 
etc., a computer simulation is an easier approach that convinced us that equation (19) is 
consistent. Although the operator A : C[0, 1] C[0, 1] 

{Ap){x) = — ---/ - r{z/x)p{z)dz (22) 

x{l — x)c{x) Jx \x J \1 — X J 

is not a contraction in the uniform topology^, the iteration procedure p„ = Apn-i con- 
verges to a smooth and strictly positive solution p{x) for all tested probability distributions 
r(s), initial approximations po, and values of e (—1, oo). Figures 2-4 show the numeri- 
cal solutions p{x) and n{s) of equation (19) for the pair-interaction R{s) = \e~'^^'^^^^ and 
"earning" P{s) = 2/(1 -|- (s — 1)'^); black, red and blue curves correspond to d = 0, d > 0, 
and d < 0, respectively. 

It is worthwhile noting that the solution for equation (19) can be calculated on any 
interval (y, 1] C (0,1]. This means that problem (6) remains well-posed on any finite 
interval [smin, Smax] and n(smax) = O^s'^l^). 

It was observed that the Fourier or Chebyshev polynomial approximation of the solu- 
tion p{x) does not provide a fast convergence. By differentiating equation (19), one can 
derive a pair of homogeneous boundary conditions of the form p'{x) — ^^pi^), x — 0,1. 
Note that the standard orthogonal polynomials do not fulfill the boundary conditions of 
this type; this is a possible explanation of the inefficiency of the polynomial approxima- 
tion. Nevertheless, the finite dimensional approximation of operator (22) in the Fourier 
basis shows that the spectrum {X^} oi I — An is distributed more or less uniformly over 
a smooth "semicircle" curve (D^eAfe e [—1,0], 3mXk G [0,1/2]) (see Fig. 5). 

^The estimate given in Remark 3 implies that \\A\\c < F,nin2 < oo, and the numerical spectral 
estimates show that \\A\\c > 1. 
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Figure 2: For ^-demand, we set a{s) = s and n{s) = no s ^(1 — s)'^p{l/s); more precisely, A = 3 



for d = (black curve) , A = 4for(i = l/3 (red curve) , A = 2 for d 
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Figure 3: In order to illustrate p and n for a-demand, we set o"(s) = 20/(1 + 19s ^). Solving 
equation (14), we find b = 2.22741 for d = and A = 3.9, b = 2.25672 for d = 1/39 and A = 4, 
6 = 2.13304 for d = -19/39 and A = 2, so that n(s) = nos'^l - p{l/ s) / a{s). 



U-demand: rho 



U-demand: n 




0.01 



0.0001 




10 20 50 100 200 



Figure 4: For u-demand, we have a{s) = l + logs and n{s) = uqs {1 — sY p{l / s) / ij{s) . Simple 
computation implies A = 3 for d = 0, A = 4 for d = 0.3, and A = 2 for d = —1/3. 
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Figure 5: The spectrums of the operator I — A„ in i-cases for d = and d = —0.5, in the 
complex plain for n = 100. There is no reason to expect the existence of a spectral gap; it is 
most likely that the limit spectrum is continuous. 

The spectrum has no spectral gap near the minimal eigenvalue Aq ~ 0, however, in 
a finite dimensional approximation at least, the minimal eigenvalue is simple. Thus, the 
spectral approximation provides no helpful hints as to the proof of the uniqueness theorem 
based on the existence of a spectral gap. In contrast. Fig. 4 shows that all eigenvalues 
A, such that |1 — A| > 1, are complex; the corresponding eigenf unctions must also be 
complex. This observation illustrates and extends the statement which will be proved in 
Theorem 3 below: a positive bounded solution of equation (6) is unique, while all other 
solutions are complex. This statement is a kind of Perron-Frobenius theorem. 

Now we must pay attention to the interesting mathematical background to the prob- 
lem (6)- (7). The consideration below is self-consistent and does not refer the reader to 
mathematical literature. 

4 Markov chains and a priori bounds for p{x) 
4.1 Transition probabilities near and 1 

Consider the Markov chain {xn} on (0, 1] with transition probability P{x, A) from x to 
A C [x, 1] given by a pair of measures describing the process in the neighborhoods of the 
points and 1: 



where a = for i- and w-demands; for s-demand, a satisfies equation (21); and Za is a 
local normalizing factor: 



Note that Zq{x) < 1 and Za{x) < 1 + I/SqRq (see (11)), but Za{x)/{1 — x)c{x) < 1 in 
any case on (0, x*]. 




(23) 




(24) 
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The associated function F{x, y) on the set x, y : 1 > y > x 



F{x,y) = < 



1 e(l) (yY if^e(^^,l], 



X c{x) r(l) \x ) ' 

is defined such that equation (19) is equivalent to the integral equation 

p{x) = E,F{x, xi)p{xi), X e (0, 1]. (25) 

As the mathematical expectation K^, x E (0, 1], is the well-defined integral operator, any 
measurable bounded solution p{x) is absolutely continuous; moreover, equation (25) can 
be rewritten as 

n 

p{x) = Ea, Jl F{Xk-uXk) p{Xn). 
k=l 

According to Lemma 1 from the Appendix, the Markov chain exits any interval (,t, y) C 
(0, 1] with probability one. Hence, if the solution p is bounded, it is continuous and 
Pi^n) p(l) = 1 as n — > oo with probability one; therefore, 

oo 

p{x) = Ea; F{Xn-l, Xn), Xq = X, (26) 
n=l 

where {xn}f is the monotone discrete-time Markov chain to be defined by (23). 

Consider a priori bounds for mathematical expectation (26). For simplicity, it is 
supposed that the probability density r{s) is uniformly bounded. Then a; < y on the 
domain of F{x, y). It is clear that F(x, y) 1, if either y ^ or x — >■ 1. Set 

= max Fix, y) 

0<x<x»<y<l ^ ' 

and consider the integer- valued exit time = n{x, x^) > from the segment (0, x^] for 
the Markov chain {xn}^lQ staring at Xq — x: < < . We have 

oo n« — 1 oo 

Yl F{Xn-i,Xn) =-F( Xn*—lj Xji^ 

) n F{Xn-i,Xn) n F{ym-l,ym) 
n=l n=l m=l 

n, — 1 oo 

- n F{Xn-l,Xn) n F{ym-l:ym): 
n=l m=l 

where {ym} is the Markov chain starting at the point y = yo = Xn» e [x*, 1]. Thus, 

oo n, — 1 oo 

Xn—lj Xn ) < F^Er, Yl F{xn-i,Xn) supEy Yl F{y^_i,y^). 

n=l n=l 3/^3;* rn=l 

Therefore, it suffices to prove the uniform boundedness of ni(a;) = ExYlnF{xn-i,Xn) for 
Xk e (0, x^] and 112(2;) = E^; Um F{yn-i, Vn) for ym e [x*, 1). 
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4.2 Upper bound and the uniqueness theorem 

Consider the uniform upper bound for IIi. According to (11) and (24), 

Za{x) / {1 — x)c{x) < 1 \/x < X^. 

The Markov chain {xn} increases monotonically and exits any interval {x,y) C (0, 
with probabihty one (see Appendix, Lemma 1). Thus, we have 



n.(.) = "n .„) < (i^^) < { „ ;[ ^ ^ °- ) 

Consider the upper bound for 112. Note that F(l, 1) = 1, and suppose that 

1-F{x,y) 

^ _ c(l) r{y/x) /y\ 
xc(x) r(l) \xj 



< oo. 



TP* def 1 

F — max 

x>x*,ye[x,l] 1 — X 
1 

max 



x>x»,yelx,l] 1 — X 



< OO. (27) 



For example, if R{s) is a monotone function and a — 0, then r{y/x) < r(l), and 

< J- fc(l) + max Mll^fM^ , (28) 
where is minc(x). In any event, F(y„_i,y„) < 1 + (1 — yn-i) F*, and the estimate 

oo 

U2{y) =^2/11 ^iym-l: Vm) < exp{F* ^(1 - y^)} , yo = y 



m=l 



can be rewritten in terms of the independent variables ^„ =^ (1 — yn)/{l — yn-i) G (0, 1] 
as the sum of a stochastic geometrical progression: 

OD / OO m \ 

J2ii-ym) = ii-y)(i+EIl^^ Ue(0,l]. 

m=0 \ m=lfe=l / 

Theorem 2 Let F* be defined by (27). Then, p{x) is a bounded continuous function if 

d+KF* F*<\n2. (29) 

Proof. As proved above, supj.g(o,s*] ni(x) < oo and < oo under rather general 
assumptions. To prove that sup^^g^^,^ n2(y) < oo, it suffices to ensure the finiteness of 
the exponential moment 

My) = Ey e(^-^)^* (^+«) < oo, G = £ g r;„. (30) 

m=l k=l 
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According to (23), the factors r]m ^ [0, 1] of the random geometrical progression (30) are 
independent random variables with the following probability distribution: 

P{r)meB)^-f [ Tj'^dr), 7 = (i+l>0. (31) 
Jb 

Note that //„ = 7/(7 + for random variable (31), and that according to (36) from 
Lemma 2 (see Appendix), the exponential moment of the geometric progression (30) is 
finite if 

supf^-^^f^fl + ^'l <e^-l + ^(e^-l-A)<l 

for X = F*. This inequality can be resolved with respect to 7: 7 < A (2 — e'^) / (e'^ — 1 — A). 
Integral equation (25) implies the absolute continuity of a bounded solution p. □ 

Remark 3 Condition (29) is the most restrictive assumption. For example, as R{s) is 
a monotone function in cases (i-iii) considered in Section 3, estimate (28) can be used to 
evaluate = 20. The convergence still holds despite the fact that the sufficient condition 

< ln2 f« 0.693 (see (29)) is crudely violated, i.e., the sufficient condition is rather far 
from an unknown necessary criteria. It appears that assumption (29) can be relaxed after 
a more detailed study of exponential moments. 

A remarkable fact to arise from this analysis is that a solution, bounded from above, 
is unique, independently of the lower bound. 

Theorem 3 The class of strictly positive bounded continuous solutions of equation (15) 
consists of the unique function. 

Proof. Suppose that there exists two nonncgative bounded continuous solutions, pi{x) 
and P2{x), that satisfy the condition p(l) = 1. The continuity assumption impUes that 
di — d2 and bi — 62 for ni{s) and ^2(5), respectively. 

Consider the bounded and uniformly continuous function A(x) = \pi{x) — P2{x)\ on 
[0, 1]. It is clear that 

max P(x,y) = — min / di = 1 — < 1 

xelx*,y] xelx»,y]Jx \1 — X J \1 — X^ J 

majorizes the probability for the Markov chain to remain in the segment [x,y], y < 1, 
X G after one step, and hence the chain exits any segment with a probability 

of one (see also Lemma 1 in Appendix). In contrast, K^H^ F{xn-i,Xn) < 112(2;) < 00, 
where 112(2;) is uniformly bounded in [a;*, 1) (see Theorem 2). Equation (25) implies that 

Uy 

A(x) < n 

Xji) max A{y) < 112 (x) max A(^), 

n=i y^l^ny,!] iely,i] 

where Uy is the exit time from [x*, y]: Xn -1 < y < Xn ■ As y can be arbitrarily chosen to 
be close to 1 and 

lim max A(^) = 0, 

(because of the uniform continuity) the previous estimate implies that A(a;) = for all 
X > x^. The uniform boundedness of Ili{y), y G (0,a;*] enables us to apply similar 
arguments to prove that A(|/) = 0, y e (0,a;*). □ 
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4.3 Lower bounds 



This subsection provides a proof for the existence of a strictly positive lower bound for 
p{x) in the i- and s- cases. For w-demand, in the particular case cr(s) = 1 + Ins, a 
weaker estimate of p{x) = 0{{liix~^)~^), g > will be established in the neighborhood 
of the origin (see Lemma 4 in Appendix). This estimate does not change the value of 
the Pareto exponent, and because of Remark 1, the necessary continuity condition 6„ = 2 
also remains valid. 

As p{x) is continuous and p(l) = 1, it follows that p is strictly positive on a certain 
segment [x^, 1], i.e., there exists such that 

min p(x) > > 0. 

a;e[x.,l] ~ 

For simplicity, suppose that the point x^, is the same as that in equations (11) and (23). 
Let the s-demand be given by a smooth increasing bounded function a{s) and that the 
tail of the probability distribution r(s) is not too heavy, i.e., there exist p such that 

Z^{X) 1 



r s°r(s) ds > — Vx e (0, x^]. (32) 



(1 — x) c{x) (1 — x) c{x) Ji 1 + p 

For i- and s-demands, this estimate holds for r(s) with any tail dominated by 0(s~°~^). 
Under the same assumption on r(s), for u-demand we have 

ZJx) If I 

^ 1 1 + D / _ix ■ (33) 



(1 — x) c{x) 1 + px \ Roa{x 1) 

In the following theorem, the lower bounds are established for p under assumptions (32)- 
(33) for any demand, but in u-c&se, as a typical demand function a{s) = 1 + Ins will 
be considered. The proof remains vahd for a wider class of slowly varying functions that 
inherit the property of the logarithm a{Y[Sk) < C Z^(7(sjk), with a constant C that is 
independent of the number of factors {sfe}. 

Theorem 4 For £- and s-demands, p{x) is strictly positive, while for u-demand, p{x) 
is bounded below by 1/ f{x~^), where f{s) a slowly varying function. More precisely, 
f{s) = (lns)» and 

p{x) >p*e-^(lnx-^)-^ 
with positive constants C and g that are specified in the proof. 

Proof. The proof does not differ for and s-cases. Lemma 1 implies the a.s. existence of 
a finite exit time n^. = n*(a;,a;*) from the segment (a;, a;*] for the Markov chain (23) {x^} 
that starts at the point Xq — x: < x^ < . Suppose that d <0. Hence, we have 



p{x) =E^'[[F{Xn-uXn) = ^x'^Xn, YlFiXn-l, X^) 

1 



n n 

n« — 1 rif — l 

n F{^n-l:Xn)Exn,-iF{Xn,-l:Xn,)p{XnJ > P* H T 

n=l n=0 



+ pXr. 
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Taking into account the inequalities Ee^ > e^^ and ln(l + x) < x, we obtain the estimate 

p{x) > expl -i^Er, J2 ^n\, X e (0, x^], (34) 

I n=0 J 

where the points Xn can be expressed as the product of factors = Xk+i/xk £ (l,oo), 
whose distribution is given by the first equation (23): 

Xn^xYl ^Xn,-1 n ( ] <X* Yl 

k=l^k-l fe=n-l\^'=-l/ k=n-l 

1 /-l/a; 



Note that the function E^.^"^ increases in a; e (0,a;*], i.e., E^.^"^ < E^,^"^ e* < 1, and 
the Markov property imphes the estimate 

n»— 1 

E^x„ = E^x„,_i n E^/^r^ <3:^* IZ< < (l-e*)~\ a;* < 1. 

fe=n-l n 

Hence, the scries E^^^ J]„ x„ converges uniformly in this fact implies the strict positivity 
of p: p{x) > p* cxp{— /i/ (1 — e*)} > 0, where ^ is defined by (32). 
Consider the case d > 0. By the Markov property, we have 

-1 X d ru-l 1 n,-l 1 /1-r 



, TT — ^ — = TT — - — A^-^) 

n* — 1 -1 

> inf E,^(l-a;^+i)'^E, [] TT ' (^5) 



where the expectation of the last product is finite (see the proof for d < 0) and the first 
infimum is strictly positive: 



inf E.Jl-Xm+i)''^ inf C'%''r{s) {I - sxf ds 



> mf ^ . . — / s r(s) ds = ^. ' , ' > 

because the integral normalized by in the last line of the equation is a decreasing 
function of x. 

Consider the it-case with (t{s) — 1 + Ins. For simplicity, suppose that d < 0; in the 
case d > the proof is similar to that in the previous paragraph. Under the assumption 
r(s) < 0(s~^), there exists // > such that 

Za{x) ^ 1 1 



(1 - x)c{x) ~ 1 + nx 1 + l/i?olnx-i 
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(see (33)). By Lemma 4, the exit time can be restricted by N[x^,x) with a probabihty 
greater than or equal to {2Y*[x)/Y^, — 1)~^, where Y*[x) = ln[x^,/x) and Y^, = E^-^ln^. 
Then, similarly to (34), we get 

p{x) >p*Ex exp<^ - iiXn + 



— i?olnx„VJ 

r N{x,x-t) 



2Y*{x)/Y^-l \ ^0 V i?olnx„V 

p* f ^('^''^*Y 1 

-2Y*{x)/Y.-l '\ V "i^olnx-^ 



where the first sum is uniformly bounded in x (see the above proof for the I- and s-cases) . 
As ln(l/x) is a concave function on x e (0, e~^], it is assumed that x^ < and apply 
the Jensen inequality to estimate the last sum: 

Note that the function E^^"^ = Sl''' r{s) ds increases in x. Hence, E-^^ ^ < 
and therefore 

111 

In — J < In h In — tj-, ^ < (N(x, x^) — n + 1) A*. 

Xn,^x\.\.k=n^k llfc=„ Ej,^^ 

where A* = ln(Ea;,^~^)~^ > 0. This estimate enables completion of the proof: 

N{x,x.) -, . N{x,x.) ^ l+hiNixx) 

i?olnx„i A^it^o ^1 n X^Rq 

where N{x, x*) < ^ In ^, i.e., the assertion of the theorem p{x) > e"*" (lnx~^)~^ holds 
with 

A* Rq 1 — A* Rq 

It is clear that the factor (Inx"-^)"^, which is a slowly varying function, does not disturb 
the Pareto exponent. □ 



5 Conclusions 

The aim of this paper is to develop a qualitative model that implies stable income distribu- 
tions with credible Pareto tails. To this end, the extreme case of economical cooperation 
has been considered and the basic features of chain-trading structures with rigid corpora- 
tive hierarchy have been axiomatized. An analysis of Section 2 shows that the qualitative 
definition of demand functions determines the values of Pareto exponents for gross and 
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net income. In general, the choice of the demand function is predetermined by the eco- 
nomic environment and the scale of the business. Hence, the qualitative behavior of the 
demand function can be motivated by estimates of market size and other macroeconomic 
constraints. The observations that are valid in terms of our model may provide some hints 
for the interpretation of facts established in empirical studies. For example, according to 
a number of reports (see [8]-[9], [17]-[18]), the graph of the net income distribution drops 
down from the Pareto exponent a = 1 to a greater value 1 + a. This phenomena can 
be explained as a transition from a long-range w-demand to s-limitation of the market 
reached by the largest corporations. 

The existence of the peak of the density distribution in the neighborhood of the mini- 
mal income {d < 0) depends on the amplitude R{1) of the pair-potential and the intensity 
of interactions Rq = J R{s) ds; the accumulation of the poverty cannot be controlled by 
the only amplitude P{1) of the social welfare. To remove the peak of poverty (to make 
d > 0) or at least reduce it, according to equation (4) it is better to decrease the gradient 
of the social aid {P'(l) 0), increase the amplitude of pair-interactions -R(l) (i.e., to 
increase the economic activity of agents), and to decrease Rq/ R{1) = J R{x)/R{l) dx (to 
decrease the range of spending). These expectable issues mean that the mathematical 
caricature (6)-(7) is potentially consistent. The quotation "// / give out my money to the 
poor in the morning, all money will be repaid by tonight" (apparently made by George 
Soros) aphoristically expresses the concept of the hierarchical model and motivates the 
basic equations of this paper. 

6 Mathematical Appendix: Four lemmas on multi- 
plicative Markov chains 

Lemma 1 (On the a.s. finiteness of the exit time) The Markov chain {xn} with transition 
probabilities (23) starting at the point Xq = x exits any interval {x,y) C (0, 1) in a finite 
number of steps a.s. 

Proof, (a) If X e (0, x*], then by definition (23), 

1 /■! 1 fl/Xn 

> (1 + £)xn\xn} = — ^7 — ^ / s'^r{y / x) dy = — — - / s"r(s) ds, 

where the right-hand side is a decreasing function of x„, as the derivative in Xn of the 
right-hand side of the equation is nonpositive. If Ji^^* r{s) ds — 0, the Markov chain exits 
{x,x^] in the first step, with a probability of one. Otherwise, set = Xn+i/xn > 1. As 
jy^* s"r(s) o?s > and r e Li, there exists £ > such that 

n^n+i >l + e}> /'^"' s°r(s) > 4 > 0, s< x^^ - 1 

independently of the previous events. The probability Pe{n, N) of the event "n random 
variables among N take values larger than 1 -|- s" can be estimated from the binomial 
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distribution: 

TV N 



P,{n, N) =p{n ^fe > (1 + er} > C% 5,^(1 - 5.) 

,4(iV + l-n) 



N-k 

k=n+l 

n 

= 1 - 0(exp{nln - A^ln(l - 5,)-^})^ 1 



as — > oo for any finite n and x < x*([24], Ch.6, §3). Hence, the Markov chain exits the 
interval {x, x*] with probabihty one. 

(b) If X e (x*, y], then by (23), for e < 2/{l + y) we have 

PK+i > (1 + e)xn\xn} = - 0? f ^— = fl + £ - T-^y 

which is clearly a decreasing function of Xn- Hence, 

PK+i > (1 + e)xn\xn} > F{x > (1 + e)y\y} ^ ^1 + e - j ^ > 

independently of the previous events. By the same arguments as in case (a), for any finite 
n and x > x*, we have Pjni^ Cfe ^ (1 + ^ 1 as A?" — > oo. □ 

Lemma 2 (On the finiteness of exponential moments) Let be i.i.r.v. with moments 
= E ^"^ < 1 and let G he the random geometric progression 

oo n 

^-Enefe=ei(i+6(i+e3(i+. ..)))■ 

n=l fc=l 



Then E e'*''-^ < oo for all X E C such that 

A: a(A) sup E ^ <1- (36) 

Proof. It sTifRccs to prove the statement in the case A > 0. Let us regard sum (30) as a 
product of two random variables, ^ and 1 + G, 

G' = 6(i + 6(i + ---)) = 6(i + G), 

where and G are independent, G and G arc identically distributed. Denote by the 
moments of the random variable G: Mn == / P{dG). Therefore, the expansion 



implies a solvable recurrent system of equations for the moments M„: 



n-l 



M„ = -^ X^C^Mfc, Mo = l, (37) 

fe=0 
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where = k\{n-ky. • Using the recurrent definition (37) of M„, we obtain 

\ fc=l / 

^ A" 1 ^ A^Mfc ~ A"-^ 1 

oo \m 



oo ji^ \ n 



n=l 



~ A'^ 1 ^ A'^ Mfe ~ A 

n=l /^n k=l m=] 



< L + 5 sup ^ 



A^ 



A;! ^1 m! 1 - 

1 



L + a{X)S, L 



k>l m=l 1 - fJ'k+m 

Thus, S <L{1- ct(A))-i < OO, provided a{X) < 1. 



\ n\ (1 - ■ 



□ 



Lemma 3 (Anti-Chcbyshev inequahty) Let {Yn} be a sequence of positive finite random 
variables such that Yn < Y* < oo and inf„^y^_...^y^ E(F„+i|Fi, . . . , Y^) = > 0. Then, 



nYlYk>7^nYA> 



1 



<k=l 



2Y*/Y,-1 



Proof. For X = ELi^fe^ we have nY* > X* > EX > nY^. Hence, for any e e (0, 1), 
we have ¥{X > enY^}> F{X > eEX}. In addition, 

EX < X*¥{X > eEX} + eEX{l-¥{X > eEX}). 

Therefore, by setting £ = 1/2, we obtain 



F{X > enY,} > F{X > eEX} > 



1-e 



> 



1-e 



X*/EX-e - Y*/Y^-e 



£=1/2 



2Y*/Y, - 1 



□ 



Lemma 4 (On the distribution of the exit time) The Markov chain {a;„} with transition 
probabilities (23) exits the interval {x, x^) in N{x, x^) steps with the probability 



F{n,{x) > N{x,x^)} > 
where Y*{x) = ln(x*/x) < Inx^^, and 



2Y*(x)/Y, - 1 



N{x, X*) 



2 , x^ 
— In — 



0{lnx-')-\ 



+ l = 0(lnx-^), Y^^ ] /^^^ lnss"r(s)ds > 0. 



Proof. Note that the function 



E^ In- 



Xr. 



1/Xn 



is a decreasing function of the initial point x„ e (0, x*]; hence, E^^ In > E^^ In^, and 
for the random variables y„ = In we have 



n = min E,{F,+i|ri, . . . , Fn} = In 
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Hence, for Xn — x Yl 



fe=i ^k/xk-i — X exp{X]fe=o l^^fe}) Lemma 3 implies 




k=0 




I fc=0 



2Y*(x)/Y,-V 



1 



i.e., exit time N{x, x*) 



A In ^ +1 has a probability of not less than {2Y* {x)/Y^ — l) ^ 

□ 
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